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Abstract 

We investigate the so-called superWIMP scenario with gravitino as the lightest supersymmetric 
particle (LSP) in the context of non-standard cosmology, in particular, brane world cosmology. As a 
candidate of the next-to-LSP (NLSP), we examine slepton and sneutrino. Brane world cosmological 
effects dramatically enhance the relic density of the slepton or sneutrino NLSP, so that the NLSP 
with mass of order 100 GeV can provide the correct abundance of gravitino dark matter through 
its decay. We find that with an appropriate five dimensional Planck mass, this scenario can be 
realized consistently with the constraints from Big Bang Nucleosynthesis (BBN) for both NLSP 
candidates of slepton and sneutrino. The BBN constraints for slepton NLSP are more stringent 
than that for sneutrino, as the result, the gravitino must be rather warm in the slepton NLSP case. 
The energy density of gravitino produced by thermal scattering is highly suppressed and negligible 
due to the brane world cosmological effects. 
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I. INTRODUCTION 



Various cosmological observations have revealed the ingredients and their abundance in 
the Universe, and showed that the standard ACDM cosmological model agrees with these 
data well. In fact, the abundance of cold dark matter is also precisely measured as 

n CDM h 2 ~ 0.11, (i) 

by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite [1]. However, the identity 
and the origin of dark matter particles are still open questions in cosmology and particle 
physics. 

The lightest supersymmetric particle (LSP) in supersymmetric models with the conserved 
R-parity is suitable for cold dark matter. Well-studied candidates are the lightest neutralino, 
gravitino, and axino. Historically the lightest sneutrino also has been considered as a dark 
matter candidate [2j], but it was excluded because of its relic abundance being too small and, 



conclusively, the null result in the direct dark matter searches js] 1 . 

Presumably, among them, the lightest neutralino is the most attractive candidate from 
various points of view. In more general terms, the weakly interacting massive particle 
(WIMP) is the attractive dark matter candidate, which the neutralino belongs to. The 
most interesting feature of the WIMP dark matter is that it can be thermal relics and its 
abundance is described as [3] 



^(moM/T^GeV- 1 



n DM h* ~ 1 x 10 ^' Tu , \ > ( 2 ) 



where mpM is the mass of the dark matter particle, Tf is the freeze out temperature, g* is 
the effective total number of relativistic degrees of freedom, Mpi = 1.2 x 10 19 GeV is the 
Planck mass, and (av) is the thermal averaged product of the annihilation cross section and 
the relative velocity. This relic abundance is essentially determined by only its annihilation 
cross section, and does not depend on the detail thermal history of the early Universe before 
T > 10 GeV. This predictability and defmiteness is fascinating, compared with the most 
of models for nonthermal production of dark matter where the resultant dark matter relic 



1 Sneutrino here means a left-handed (active) one in the minimal supersymmetric standard model (MSSM) 
with the lepton number conservation. The case with the lepton number violation is discussed in Ref. 
A right-handed (sterile) sneutrino [5[ or a mixture of left- and right-handed sneutrino [6| is still viable. 



2 



density is highly model dependent. In addition, on-going or forthcoming experiments may 
observe such a dark matter directly and/or indirectly. 

However, there is another interesting possibility for the candidate of dark matter with 
highly suppressed interactions. Even though it cannot be in thermal equilibrium, this type of 
dark matter can have the appealing feature in its abundance. This is called the superWIMP 
scenario, where the superWIMP dark matter is produced through the decay of a long-lived 
WIMP [8|, |9[ with highly suppressed interactions between the WIMP and the superWIMP. 
In this case, the dark matter abundance is given as 

n DM h 2 = ^n x h 2 , (3) 

m x 

where Qxh 2 would be the thermal relic abundance of the WIMP given by Eq. if it were 
stable. 

There are two typical superWIMP scenarios in the context of supersymmetric models. 
One is a scenario realizing the axino as the LSP with interactions suppressed by the Peccei- 
Quinn (PQ) symmetry breaking scale 8| . The other is a scenario where the gravitino is the 
LSP with Planck scale suppressed interactions 

The superWIMP scenario is an interesting possibility; however, it is not so easy to consis- 
tently realize the scenario with big bang nucleosynthesis (BBN). In supersymmetric models, 
the next-to-LSP (NLSP) of, e.g., neutralino, stau, sneutrino decays into the LSP and the 
standard model (SM) particles at late times. If the NLSP decays after BBN, its energetic 
daughters of the SM particles would destroy light nuclei and spoil the success of the standard 
BBN predictions. Furthermore, late-time injection of the energetic photon from the NLSP 
decay would distort the spectrum of the observed cosmic microwave background. 

In the case of the axino LSP dark matter, the constraints from BBN can be rather easily 

n 

avoided because the NLSP has a relatively short lifetime [8J and decays well before BBN. 
Some variants of this scenario are also possible [lo| . 

In contrast to the axino LSP scenario, for the case of gravitino LSP dark matter, the 
constraints from BBN turn out to be very stringent. This is because the NLSP is normally 
long-lived through the Planck mass suppressed interaction between the NLSP and the grav- 
itino LSP. When the NLSP is the neutralino or the stau (or the slepton in general), this 

" n n 



scenario is very hard [ill, Haiia, H4 

In the stau NLSP case, so-cal 



the bound state with charged NLSP makes the problem worse 
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ed catalyzed-BBN due to 



131 ] . On the other hand, 



the sneutrino NLSP has rather mild constraints from BBN compared with the others. Nev- 
ertheless, the sneutrino NLSP has not been appealing from the beginning in this framework 
because the thermal relic abundance Q„h 2 is, as is well-known, too small to be suitable in 
Eq. due to the very efficient annihilation, unless the NLSP sneutrino mass is of order 
of TeV j^]. Such a heavy NLSP is not appealing from the weak scale supersymmetry point 
of view. In order to obtain a suitable relic density of sneutrino NLSP in E q. (13]) , one may 
consider a nonthermal production of sneutrino NLSP (see, for example, Ref. |15|). However, 
in such a nonthermal scenario, Qxh 2 generally depends on many parameters in the scenario. 
Here, we stress that the sneutrino NLSP scenario is not suitable as the super WIMP scenario, 
because of its too small thermal relic abundance with a mass of (9(100 GeV), rather than 
the BBN constraints. 

In general, if we abandon providing the gravitino LSP relic density from the NLSP decay 
through the formula in Eq. ([3]), we may be able to avoid the BBN constraints in some 
parameter region in a model. In this case, the gravitino LSP should be produced through 
scattering and decay processes of particles in thermal plasma (thermal pro duction). The 



relic density of the gravitino through thermal production is evaluated as 16 1 



where Tr is the reheating temperature after inflation, and M 3 is the running gluino mass. By 
adjusting the reheating temperature, gravitino mass, etc., we can obtain the relic abundance 
consistent with the WMAP data. However, this scenario may not be appealing because the 
dark matter relic density has nothing to do with the NLSP relic density. 

In this paper, we consider the superWIMP scenario with the gravitino LSP in the context 
of nonstandard cosmological modes with a modified expansion rate of the early Universe. 
The nonstandard cosmological effects cause dramatic changes into the thermal relic density 
of the NLSP with mass of (9(100 GeV), so that this NLSP decay can yield the correct 
amount of the gravitino dark matter relic density. We find that this scenario can be real- 
ized consistently with the successful predictions of BBN. To be concrete, we consider the 



brane world cosmology [12] based on the picture proposed by Randall- Sundrum [18|. Other 



cosmological models such as quintessence scenarios 19] and scalar-tensor theories 20| have 
similar effects. In the brane world cosmology, we examine two possibilities for the NLSP, 
stau (or slepton, in general), and sneutrino. 



II. BRANE WORLD COSMOLOGY 



Here we give a very brief review on our nonstandard cosmological model, namely, the 
brane world cosmology, which has been intensively investigated [17[. This model is a cos- 
mological version of the so-called RS II model proposed by Randall and Sundrum [18|. In 
the model, our four-dimensional universe is realized on the 3-brane with a positive tension 
located at the ultraviolet (UV) boundary of a five- dimensional anti-de sitter spacetime. In 
this setup, the Friedmann equation for a spatially flat spacetime is given by 

3 V Po, 

where 



(5) 



po = 96nGM, 



5 ! 



(6) 



with M 5 being the five- dimensional Planck mass. Here the four-dimensional cosmological 
constant has been tuned to be zero, and the negligible so-called dark radiation has been 
omitted 2l|. The second term proportional to p 2 is a new ingredient in the brane world 



cosmology and leads to a nonstandard expansion law of the early Universe. In the following, 
we call a temperature T t defined as p{T t ) = |^-g*T 4 4 = p "transition temperature." This 
modification of the expansion law H 2 oc p 2 at a high temperature (T > T t ) leads to some 
drastic changes for several cosmological issues, not only the thermal relic abundance of dark 



matter [2J 



23 



24 



251 ] but also other subjects, e.g., leptogenesis [26[]. Needless to say, 
until the onset of BBN, the standard expansion law, H 2 oc p, has to be recovered. This 
requirement leads to the model-independent constraint on the transition temperature as 
T t > 1 MeV or equivalently M 5 > 8.8 TeV [17|. 

To be precise, we specify our setup here. We discuss a super symmetric dark matter 
scenario, and our model should be a supersymmetric version of the Randall-Sundrum brane 

n 

world model [27j. Since the Einstein equation belongs to the bosonic part in supergravity, 
the cosmological solution of Eq. (jHJ) remains the same. Although the gravitino is a bulk 



field, we can, as a good a 
in the previous works 28 . 



jproximation, regard it as a field residing on the UV brane as 



29 



30t ] because the zero-mode gravitino is localized around the 



UV brane as the same as the zero-mode graviton is. We assume the cancellation of the 
four-dimensional cosmological constant by some mechanism different from that discussed 



5 



in the original paper [18j, where the model parameter is severely constrained by precision 
measurements of the gravitational law in sub millimeter range, T t > 1.3 TeV or equivalently 
M 5 > 1.1 x 10 8 GeV. This result is quite model dependent. For example, in dilatonic brane 



world models 31], which seem to be a more natural five-dimensional effective theory of an 
underlying fundamental theory than the original RS model, the constraint can be very much 
moderated depending on the stabilization mechanism of the dilaton field. Therefore, we take 
only the BBN bound into account to keep our discussion general. Moreover, in this paper, 
we do not specify the origin of supersymmetry breaking and its mediation to the visible 
sector, so that we treat the NLSP mass of 0(100 GeV) and the gravitino mass as a free 
parameter. 

III. GRAVITINO ABUNDANCE FROM ENHANCED RELIC DENSITY 



As has been pointed out in Ref . [22| , the relic density of dark matter can be enhanced in 
brane world cosmology when the five-dimensional Planck mass M5 is low enough. Applying 
this effect to the abundance of NLSPs, we calculate the relic density of the gravitino LSP 
produced by the decay of the NLSP at late times. 

In the context of brane world cosmology, the thermal relic density of a stable or a long- 
lived particle by solving the Boltzmann equation, 

^ + 3Hn = -(av)(n 2 -nl q ), (7) 

with the modified Friedmann equation, Eq. (jSJ), where n is the actual number density of the 
particle, and n eq is the equilibrium number density. If the particle freezes out at the brane 

we can solve the Boltzmann equation 



world cosmology era where p 2 dominates in Eq. J5 
analytically and find an approximation formula 22 

Y = - * 0.54-f^. (8) 
s \{av) 

1 /9 

Here A = 0.26(g* s /g*' )M Pl m, m is the mass of the particle, xt is defined as x% = m/T t , 
and we have assumed the annihilation processes in S-wave. Note that the relic abundance 
is characterized by the transition temperature rather than the decoupling temperature. 
Using the well-known approximate formula in the standard cosmology 



with the freeze-out temperature Xd = fr, we obtain the ratio of the relic energy density of 
dark matter in the brane world cosmology to the one in the standard cosmology {Q( s )) 

such that 

^ - 0.54 (*_) , ( 10 ) 

where Xd( s ) denotes the freeze-out temperature in the standard cosmology. Therefore, the 
relic energy density in the brane world cosmology can be enhanced from the one in the 
standard cosmology when the transition temperature is low enough, T t < Tf. In the following 
analysis, we assume such a low transition temperature. 



A. Slepton (stau) NLSP 

First, let us consider the case of the slepton NLSP. Among many possible channels of 
annihilation processes, the dominant one is the annihilation into 77 and Z'j. Hence, the 



cross section can be roughly given as 



32] 



47ra 2 



With this annihilation cross section, the re 



the standard cosmology is expressed as 



32] 



ic abundance of the slepton (if it were stable) in 



m,-yr-10- n f m X VGeV, (12) 
' ' VIOOGeV/ ' V ' 

or in terms of the density parameter, 

Using Eqs. fflUj) . ffT2|) and (TT3T) . the slepton NLSP abundance in the brane world cosmology 
is modified into 

mt Y m = 4 x KTlGeV] x (^J' (JL) (3.) , (M) 



Thus, the gravitino abundance from the slepton LSP decay is given as 

' 23 ' 

X d(s), 



3/2 mj m \ mj J VIOOGeV/ \x d(a) ) V1700/ V ; 



Although appropriate parameters enable us to obtain the correct relic density of the dark 
matter gravitino, these parameters are severely constrained by BBN. The consistency with 
the constraints from BBN is the key whether the superWIMP scenario with the gravitino 
LSP can be realized or not. 

Now let us consider the constraints from BBN. Detailed analysis for the stau NLSP case 
has been done, and the results are summarized in Figs. 2 and 3 for m^„ = 100 GeV and 300 

n 

GeV, respectively, in Ref. [13J] . In fact, we can see a large parameter region is excluded. For 
nif R = 100 GeV, the region, which is consistent with BBN and provides a suitable amount 
of the gravitino relic density via the stau NLSP decay, can be read off as 

(m fR Y fR ,m 3/2 ) ~ (10- 6 GeV, 10 MeV). (17) 

In order to realize this allowed region, Eq. (fl4|) indicates T t ~ 24 keV for xm s ) ~ 23, so 
that this result contradicts against the requirement T f > 1 MeV. For m^- R = 300 GeV, the 
allowed region appears in 

m ffl y fB >2xlO- 7 GeV, m 3/2 <0.6GeV. (18) 

We find that the transition temperature, T t < 3.2 MeV (or equivalently M 5 < 19 TeV), can 
realize this region. This result is close to the BBN bound on the transition temperature 
T t > 1 MeV, and the gravitino superWIMP scenario with the stau NLSP is marginally 
possible in the brane world scenario. Note that the standard cosmology cannot realize this 
scenario as can be seen from Eq. ({TBI) and the enhancement of the stau NLSP abundance 
by the brane world cosmological effects is necessary. 

In the allowed region we have considered, the gravitino is as light as sub-GeV while the 
NLSP is as heavy as 100 GeV. Hence, this gravitino dark matter produced by the decay of 
stau may have a large velocity dispersion. While cold dark matter (CDM) has made great 
success in the structure formation of the Universe at a large scale (> 1 Mpc), the recent high- 
resolution iV-body simulations of the CDM showed some discrepancies at a small scale (< 1 



Mpc). One is the so-called "missing satellite problem" [33j; namely, the iV-body simulation 
predicts more virialized dark objects with much smaller mass scale than observed. The 
other is that CDM models predict a cuspier mass distribution of CDM halos than observed: 
"cusp problem" [3^. It has been discussed 35, 36] that these problems can be solved in the 



superWIMP scenario with appropriate large velocity dispersion of the superWIMP, which 



reduces the number of less massive subhalos and the primordial fluctuations at the small 



scales; its difficulty is also pointed out though [37]] . On the other hand, too large velocity 
dispersion of the superWIMP dark matter, which exceeds the damping scale > 1 Mpc, is 
constrained from Lyman alpha clouds. The damping scale can be calculated for a given 



lifetime of NLSP 



38] 



2 

T-NLSP ^ 48ttM£-^-, (19) 
m NLSP 

with Mp being the reduced Planck mass, the NLSP mass m m&y and the gravitino mass 
m 3 / 2 , and a given mass ratio m 3 / 2 /^NLSP ( see Fig- 1 i n Ref. |36|). In our case with m 3 / 2 
in sub-GeV and slepton NLSP mass around 100 GeV, typical damping scale is found to be 
(9(1 Mpc). Hence, our scenario is in fact marginally possible. 

In our analysis, we have only considered the gravitino abundance from the decay of 
the NLSP. In general, the gravitino production from thermal plasma should be considered. 



As has been pointed out in Ref. 



281 ]. the relic density of the gravitino through the thermal 



production in the brane world cosmology can be obtained from Eq. with the replacement 
Tr — > 2T t in the case with T t < Tr, 

— (SG3(^) 2 

The gravitino density from the thermal production is negligible unless Tr > 10 8 GeV. This 
is also an important effect in the brane world cosmology. The gravitino abundance given in 
Eq. (fl6|) is reduced into the standard cosmological one for x t < Xd( s ), 



^O^O-^j^j, (21) 

where we have used Eq. f fl3|) . For m z ~ ~ 100 GeV in our interest, the gravitino relic abundance 
in the standard cosmology is too small to be consistent with the WMAP data as mentioned 
above. As the transition temperature is raised from T t ~ m ; -/1700 [see Eq. ffT6l) ]. the total 
gravitino relic density from the sum of the NLSP decay and the thermal production is 
decreasing and reaches the standard cosmological value. When the transition temperature 
is raised further, the gravitino relic density from the thermal production exceeds that from 
the NLSP decay for T t > 10 6 GeV and becomes consistent with the WMAP data at T t ~ 10 8 
GeV. 



2 
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B. Sneutrino NLSP 



n the standard cosmology, the density parameter for the relic sneutrino is estimated 
a s [2I 

O^.O^xf^) 2 . (22) 
Together with Eqs. (fTUj) and ([221 . we find, in the brane world cosmology, 

n m tf * 0.1 x ( m > V ( (-*-) , (23) 
( ) VlOOGeV/ \x d{s) J V2100/ ' 1 ; 

or equivalently 

Therefore, due to the enhancement by the modified expansion law, the sneutrino with mass 
0(100 GeV) can provide a suitable relic density of the gravitino LSP via its decay. The 
resultant density parameter of gravitino is given as 

n 3/2 h* = * 0.1 x ( ^) ( m ; ) 2 ( (25) 

61 \m^J\100GeVJ \x d{s) J V2100/ V ; 

from Eq. (j3J). 

The BBN constraints for the sneutrino LSP have been analyzed in detail, and the results 



are summarized in Fig. 5 for niy = 300 GeV in Ref. 14j . The allowed region can be read 
off as 

(i) mpY m > 2 x 10- 6 GeV, m 3/2 < 0.06 GeV, 

(ii) m„Y m ~2x 10~ 7 GeV, m 3/2 ~ 0.3 GeV, 

(iii) 4 x 10- 10 GeV < m P Y m < 7 x lO" 10 GeV, 300 GeV > m 3/2 > 150 GeV. 

For each allowed region, Eq. (1211) leads to 

(i) T t < 0.26 MeV (M 5 < 3.6 TeV), 

(ii) T t ~ 2.6 MeV (M 5 ~ 17 TeV), 

(iii) 1.3 GeV > T t > 730 MeV (1000 TeV > M 5 > 710 TeV). 

Region (i) is not consistent with T t > 1 MeV and hence excluded. Point (ii) is marginally 
acceptable as in the case of the stau NLSP but similarly delicate from viewpoints of both 

10 



BBN and free-streaming constraints. Finally, region (iii) is typical for the sneutrino NLSP 
case, because this appears when masses of the sneutrino NLSP and the gravitino LSP are 
close and such a region is never allowed for other NLSP cases. In addition, one may notice 
that the free-streaming length of region (iii) is less than a Mpc scale unlike point (ii). 



IV. CONCLUSION AND DISCUSSIONS 



We have investigated that gravitino dark matter produced by the decay of thermal relic 
NLSP, in particular, stau NLSP and sneutrino NLSP, in the context of brane world cosmol- 
ogy 2 . Even if the NLSP mass is around 100 GeV, the brane world cosmological effects can 
dramatically enhance thermal relic density of the NLSP enough to yield the correct abun- 
dance for the gravitino dark matter through the NLSP decay. We have found that with a 
suitable five-dimensional Planck mass, this scenario can be consistent with BBN constraints. 
The brane world cosmological effects suppress the thermal production of the gravitino, so 
that the energy density from the thermal production is negligible in our scenario. 

In the brane world cosmology, a disfavored candidate of dark matter due to too efficient 
annihilation in the standard cosmology can be a suitable one with the help of the enhance- 
ment effect for its relic density. The case of the winolike neutralino dark matter has been 



analyzed in detail in Ref . 24j , and it has been found that the winolike neutralino with mass 
around 100 GeV can be consistent with the WMAP data. We have shown that the same 
enhancement mechanism is applicable to the NLSP in the superWIMP scenario with the 
gravitino LSP. The sneutrino has no electric charge and is a potential candidate of the dark 
matter, if its relic density is enhanced as in the case of winolike dark matter in Ref. 24]. 
However, the possibility of sneutrino dark matter has been excluded by the null result in 
direct dark matter searches. The sneutrino NLSP in the superWIMP scenario would be 
interesting for the following reason: There is no relic sneutrino in the present Universe and 
the sneutrino cannot be observed in direct dark matter searches, nevertheless the gravitino 
dark matter inherited its relic density from that of the sneutrino. In addition, the sneutrino 
is a harmless mother particle decaying into gravitino dark matter rather than others like the 



2 Recently, the gravitino dark matter in brane world cosmology have been investigated in Ref. [30| . In the 
paper, the NLSP is identified as the neutralino in the framework of the constrained minimal supersym- 
metric standard model (CMSSM) and the thermal production of gravitino is also taken into account. 
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neutralino or slepton, because of its weak constraints from BBN. 

When the transition temperature or equivalently M5 is very low, the enhancement of 
thermal relic density of a particle would be much too to be consistent with the WMAP 
data. This fact implies a lower bound on T t or M 5 , which can be stronger than the model- 
independent BBN bound, T t > 1 MeV or equivalently M 5 > 8.8 TeV. In fact, the lower 
bound on M5 > 600 TeV has been reported for neutralino dark matter in the CMSSM 23), 
while M5 > 100 TeV for the winolike dark matter [24J] . In the superWIMP scenario we have 
discussed, the factor, m 3 /2/m NLS p, in Eq. ([3]) can play an important role so as to moderate 
too much enhanced NLSP energy density and as a result, the lower bound on T t or M5 can 
be as low as the model-independent BBN bound. Of course, if the factor m^^/m^isp is 
too small, the resultant dark matter becomes too warm as we have discussed. 

The long-lived NLSP is also interesting at future collider experiments such as the Large 
Hadron Collider (LHC) and the International Linear Collider (ILC). In the superWIMP 
scenario with the gravitino LSP, the decay length of the NLSP well exceeds the detector size 
of the LHC and the ILC, and the NLSP decay takes place outside the detector. In this case, 
there have been interesting proposals [391] for the way to trap quasistable NLSPs outside the 
detector, when the NLSP is a charged particle like stau. Detailed studies of the NLSP decay 
may provide precise measurements of the gravitino mass and the four-dimensional Planck 
mass. In our brane world scenario, such precise measurements can also give the information 
of M 5 . 

If the long-lived NLSP is a charge neutral particle such as the neutralino and the sneu- 
trino, it is very hard to observe its decay. However, in this case, it is important to distinguish 
the sneutrino from the neutralino. At collider experiments, the long-lived sneutrino behaves 
similarly to the LSP (or long-lived NLSP) neutralino, but there are some characteristic sig- 
natures [40]. Interestingly, once the sneutrino has been observed as if it were the LSP at 
collider experiments, the superWIMP scenario would become more plausible as the dark 
matter scenario. This is because a simple scenario with the sneutrino LSP has already been 
excluded by direct dark matter search experiments and the sneutrino must be unstable, if 
the LSP is the dominant component of dark matter. As in the case of stau NLSP, precise 
measurements of sneutrino properties can be interpreted into the information of M 5 . 

If the neutralino has been observed as the LSP at collider experiments and its properties 
have been precisely measured, one can calculate the neutralino relic density at the present 
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Universe. Suppose that the result shows too large relic density inconsistent with the WMAP 
data. In this case, the superWIMP scenario can provide a solution to this contradiction: 
The neutralino is the NLSP and decays into the gravitino LSP, whose energy density is 
reduced by the factor, m 3 /2/m NLS p, in Eq. fl3j). However, the key to realize this scenario 
lies on how to avoid the constraints from BBN. 
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